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ABSTRACT 

Context. Ultra-cool dwarfs of the L spectral type (T^ff = 1400-2200 K) are known to have dusty atmospheres. Asymmetries of the 
dwarf surface may arise from rotationally-induced flattening and dust-cloud coverage, and may result in non-zero linear polarisation 
through dust scattering. 

Aims. We aim to study the heterogeneity of ultra-cool dwarfs' atmospheres and the grain-size effects on the polarisation degree in a 
sample of nine late M, L and early T dwarfs. 

Methods. We obtain linear polarimetric imaging measurements using FORS 1 at the Very Large Telescope, in the Bessel / filter, and 
for a subset in the Bessel R and the Gunn z filters. 

Results. We measure a polarisation degree of (0.31 ± 0.06)% for LHS102BC. We fail to detect linear polarisation in the rest of our 
sample, with upper-limits on the polarisation degree of each object of 0.09% to 0.76% (95% of confidence level), depending on the 
targets and the bands. For those targets we do not find evidence of large-scale cloud horizontal structure in our data. Together with 
previous surveys, our results set the fraction of ultra-cool dwarfs with detected linear polarisation to 30^^"% (l-cr errors). From the 
whole sample of well-measured objects with errors smaller than 0. 1%, the fraction of ultra-cool dwarfs with polarisation degree larger 
than 0.3% is smaller than 16% (95% confidence level). 

Conclusions. For three brown dwarfs, our observations indicate polarisation degrees different (at the 3-cr level) than previously 
reported, giving hints of possible variations. Our results fail to correlate with the current model predictions for ultra-cool dwarf 
polarisation for a flattening-induced polarisation, or with the variability studies for a polarisation induced by an heterogeneous cloud 
cover This stresses the intricacy of each of those tasks, but may as well proceed from complex and dynamic atmospheric processes. 
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1. Introduction tation induces a flattening at the poles. Furthermore, the fast 



rotation may increase the turbulence in their convective at- 
The L spectral type covers effective temperatures of 2200 K ^ospheres, leading to possible heterogeneities in the cloud 
to ^ 1400K (Ki rkpatricketalJ | 1999|: [ Martin et al. || 1999t ). The ^^^y^ ^he polarisation arising frommany small-scale, randomly- 
lowest-mass stars on the main sequence as weU as substellar distributed structures would cancel out, but large-scale hetero- 
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H . objects with a large range of masses, at a given time of their ggneities may result in net polarisation. Gravitational pulling by 

. . ^ , coohng sequence, populate this spectral type. The optical spec- ^ ^j^gg companion could also break the symmetry of the brown 

tra of L dwarfs show specific features that define the type: strong ^warf surface, as it is expected rotational acceleration does, 

hybrid bands (FeH and CrH) and alkali lines (Na i, K i, Rb i, , . ■ . 

Cs I). The oxides TiO and VO produce weak absorption fea- I Sengupta & Krishan| (12001D calculated the polarisation of 

tures, which is interpreted to be due to the coagulation of those rotationally-flattened ultra-cool dwarfs and predicted that they 

species into dust grains. Other refractory elements such as Al, '^o^l^ have non-zero polarised light. They estimated the po- 

Ca, Fe, and Mg condensate into grains and form dust clouds. By lansation degree due to flie rotationaUy-induced flattening for 

contrast, in the warmer atmosphere of the early M dwarfs, the single scattering and various grain-size distributions, while 

temperature is too high for the molecules of e.g. TiO or VO to Sengupta & Kwok (2005) studied the effect of multiple scatter- 

aggregate into grains, while the atmospheric models show that ^'^8- 

the clouds made of larger dust grains sediment below the photo- Later, two observational studies took on to test those 

sphere of the cooler mid-T brown dwaifs. predictions. Menard etal. (2002) (M02 hereafter) used the 

Several studi es have showri that ul tra-cool dwa r fs are VLT/FORSl imager in Paranal Observatory to observe in the 

fast rotating tMohantv & Basril 120031: iBailer-JonesI l2004t Bessel / band a sample of eight LO to L8 ulti-a-cool dwarfs 



IZapatero Osorio et al.ll2006h . It is therefore likely that the ro- (and binaries), detecting a significant, although small polai"- 

isation degree in three of those, between 0.1% and 0.2%. 

Send offprint requests to: B. Goldman, goldmanOmpia.de |Zapatero Osorio et al. (2005) (ZO05) used the 2.2m/CAFOS im- 

* Based on observations collected at the European Observatory, ager in Calar Alto Observatory to observe a larger sample of 

Paranal, Chile, under programmes 075.C-0653(A) and 077.C-0819(A). 43 ultra-cool dwarfs with spectral types between M6 and L7.5, 
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in Johnson R and / bands. They found eleven dwarfs with signifi- 
cant polarisation, including seven with an estimated polarisation 
degree larger than 0.5%. 

Given its dependency on the rotationally-induced flattening, 
the grain-size distribution and the large-scale cloud distribution 
on the dwarf surface, polarimetric measurements have the poten- 
tial to shed light on many fundamental processes in ultra-cool 
atmospheres, which are currently hardly constrained by other 
observing modes, or not at all. At the same time, the entangled 
dependency of the polarisation degree on those many parameters 
requires a specific observing strategy. Rotational velocities could 
be correlated to the polarisation degree. Dynamical weather-like 
cloud patterns could be detected through polarimetric monitor- 
ing. Finally, the dependence of the polarisation degree to the 
grain size can be probed through multi-wavelength observations. 
Comparison of optical and near-infrared data, particularly if the 
typical grain size is in the 1-yum range, is promising. 

This program constitutes an ambitious and demanding ob- 
serving challenge. In this article we report on the first, likely 
easiest stage, namely polarimetric observations in the red part of 
the spectrum, in one to three filters depending on the targets. 
As some targets were previously observed by M02 or ZOOS, 
multiple-epoch observations are already available. In Section|2l 
we describe our sample and the observations, as well as our data 
reduction and polarisation determination. We present our results 
in Section[3] and we compare and combine them with previous 
measurements. Finally we give our conclusions in Section|5] 

2. Observations 

2.1. Target selection 

We searched the compendium DwarfArchive^ and chose our 
targets according to the following criteria: 

1 . Availability of published polarimetric measurements, in or- 
der to search for long-term variability. 

2. Spectral type: we selected targets sampling a large range of 
spectral types from M9 to early T. 

3. Brightness: in order to achieve signal-to-noise ratios (SNR) 
of 500 and more in each image in a few minute exposure 
times, in a wide range of observing conditions, we selected 
the brightest targets in each spectral type. 

4. If possible, availability of high-resolution spectroscopic de- 
termination of the rotation velocity. 

The targets are described in TablefllFl 

2.2. Observing conditions and instrumental set-up 

All data were obtained using the Focal Reducer and l ow dis- 
persion Spectrograph FORSl (lAppenzeller et al.lll998l) . which 
is mounted at the Cassegrain focus of the UT2 (Kueyen) of 
the VLT facility. With some exceptions, we used the 2k x 2k 
Tektronix CCD and four output amplifiers. In some cases the 
CCD was windowed to the central Ik x Ik area. In both cases no 
binning was applied and the pixel scale is 0'.'2, and the full field 
of view is 6.8 x 6.8 min^. The ipol polarimetric imaging mode 
of FORSl uses a Wollaston prism splitting the beam into the 

' http : //www . Dwarf Archives . org | 

^ We used the Ultracool Dwarf Catalog hosted at 
|http : //www . iac. es/galeria/ ege/catalogo_espectral/ to 
obtain the / magnitude or estimate it from ultra-cool dwarfs of similar 
spectral type. 



ordinary and extraordinary parts; a rotating half- wave retarder 
plate; and a 11 "-wide strip mask to prevent overlapping of the 
two beams. All the targets were centred at the same CCD posi- 
tion, in the middle of the central strip. Observations were carried 
out through the Bessel R and / filters, as well as the Gunn z fil- 
ter, with central wavelength and full-width at half maximum of 
657 + 150 nm, 768 + 138 nm, and 910 + 131 nm respectively. 

We used four retarder angles: Odeg, 22.5 deg, 45deg, and 
67.5 deg, which allow to derive the Stokes parameters (see 
Section |273T l. All science targets were observed several consec- 
utive times at each retarder angle, while the bright polarimetric 
standard stars were observed only once (see Table|2|. 

The observations were conducted in two programmes, one 
in visitor mode, the second in service mode (see Table|2]i. The 
atmospheric transparency was good at all times. The seeing con- 
ditions were average for the site during service mode, but they 
were worse than average during the two visitor nights, with a 
median seeing of 1.5", and strong winds from the South-East. 

2.3. Polarimetric measurement and error determination 

The images were bias-subtracted and flat-fielded with twilight 
flat-field of the same or contiguous nights, obtained witho ut the 
Wollaston prism, using the standard noao/iraf routines dTodvl 
1993). Cosmic rays and isolated bad pixels were removed from 
the flat-field frames using the crreject iraf routine, and from 
the science images using the sigma.filter idl routine. The 
flat-fielding is an important issue, as the flat-fielding inaccura- 
cies (at the extraordinary position relative to the ordinary po- 
sition) will propagate linearly, to first order, into the Q and U 
measurements. This problem is mitigated in our case because all 



Table 2: Observing log. 



Target 


Date 


Filt. 


Seeing' 
["] 


SNR2 


^exp 

[s] 




DENIS J2200 


01/08/05 


/ 


1.7 


720 


160 


3 


Kelu 1 AB 


02/08/05 


/ 


1.4 


529 


200 


2 




02/08/05 


z 


1.6 


243 


200 


2 


2MASS J0036 


01/08/05 


R 


1.5 


706 


600 


3 




01/08/05 


I 


1.7 


677 


150 


3 




10/08/06 


I 


1.3 


785 


150 


4 




11/08/06 


I 


0.9 


937 


150 


4 


2MASS J2158 


01/08/05 


I 


1.2 


189 


300 


3 




02/08/05 


z 


1.0 


350 


400 


3 


LHS 102BC 


01/08/05 


/ 


1.2 


186 


130 


2 




02/08/05 


/ 


0.9 


192 


180 


2 




02/08/05 


z 


0.9 


532 


200 


3 




28/07/06 


I 


1.2 


890 


150 


4 


2MASS J 1507 


01/08/05 


R 


1.0 


491 


200 


3 




01/08/05 


I 


1.0 


783 


160 


3 


DENIS J2252 


02/08/05 


I 


1.0 


428 


600 


3 


2MASS J2244 


02/08/05 


z 


1.1 


124 


600 


2 




19/08/06 


/ 


0.8 


222 


300 


3 




21/08/06 


/ 


0.9 


214 


300 


3 




22/08/06 


/ 


0.8 


229 


300 


3 


£ Indi Bab 


01/08/05 


/ 


1.2 


422 


150 


3 



' The seeing is calculated from the average FWHM of the science 
target in the given night. ^ Signal-to-noise ratio (SNR), averaged for 
each target over a night. ^ Number (Ng) of exposures per retarder 
angle. 
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Table 1: Summary of targets and references, ordered by spectral type. 



Object 


1 


J 




Sp.T. 


Discovery paper 


V sin / [km/s] 


DENIS J22000201-3038327 


16.7 


13.44 ± 


0.03 


M9+L0 


Kendall et al. (2004) 


17' 


Kelu 1 AB 


16.94 


13.41 ± 


0.03 


L2 


Ruiz etal. (1997) 


60^ , 76' 


2MASS J00361617+1821104 


16.11 


12.47 ± 


0.03 


L3.5 


Reid et al. (2000) 


36±2.7\45' 


2MASS J21580457-1550098 


18.5' 


15.04 ± 


0.04 


L4 






LHS 102BC (GJ lOOlBC) 


16.68 


13.11 ± 


0.02 


L4.5+L4.5'' 


Goldman et al. (1999) 


42' 


2MASS J15074769-1627386 


16.65 


12.83 ± 


0.03 


L5 


Reid et al. (2000) 


27.2^ 32' 


DENIS-P J225210.7-173013 


17.9 


14.31 ± 


0.03 


L6:+T2:« 


Kendall et al. (2004) 




2MASS J22443 1.67+204343.3 


20.4' 


16.48 ± 


0.14 


L7.5 


Dahn et al. (2002) 




6 Indi B 


16.6 


11.91 +0.02 


T1+T6 


^cholz et al. (2003) 





' jReiners & Basri (2008): ^ Mohantv & Basril a2003i) : ^ iZapatero Osorio et alj ( 12006) ; iGolimowski et al.i ( l2004h : 
' iBailer- Jones ( 2004) : ^ Reid et al. (200^ : ^ Estimates based on ultra-cool dwarfs of similar spectral type. 



frames of the targets and standard stars are observed at the same 
CCD pixels, at an accuracy better than the full width at half max- 
imum (FWHM) of the stellar profile. The effects of the flat-field 
photonoise and systematic errors on each pixel's sensitivity are 
reduced by the excellent sampling of the point-spread function 
(PSF). A large number of pixels (typically several dozens) are 
illuminated by the star and used to measure its flux, so that the 
non-correlated effect roughly scales down as the square-root of 
that number. 

Correcting for the linear polarisation of usual flat-field 
sources including twilight sky and internal screens is difficult, 
so we obtained our flat-field frames without the Wollaston prism. 
Therefore they do not account for the whole light path. In fact, 
we measure a ^ 0.7% systematic difference between the ordi- 
nary and extraordinary fluxes in all our observations, indepen- 
dent of the filter, observing conditions, integration time or target 
magnitude. Regardless of the cause of this systematic difference, 
it mostly vanishes as for the polarisation determination we take 
the differences of the relative differences of the ordinary and ex- 
traordinary fluxes (see Eqs. below). In addition, the ordinary and 
extraordinary rays swap positions for certai n pairs of rotator po- 
sitions (when the angles differ by 7r/4, see iPatat & Romanieliol 
[2006). Ultimately the effects of the inaccuracies of the flat-field 
are estimated using the polarimetric standard stars' observations 
(see next subsection). However the spectral types and spectral 
energy distributions of the scientific targets and the standard 
stars greatly differ. 

We derive the polarimetric information in a way similar to 
that of ZOOS. We extract 30 x 20-arcsec^ stamp images around 
the target, with some adjustments when a nearby star could af- 
fect the background measurement. The fluxes of the ordinary and 
extraordinary stamp images are measured w ith aperture photom- 
etry using using SExtractor dBertin & Arnouts 1996) . as well 
as PSF fitted photometry using the iraf daophot program. For 
consistency, we chose a constant aperture of lOpixels, or 2.0". 
We then average the (extra)ordinary fluxes of the 2-A images at 
a given retarder angle, and use the standard deviation as a proxy 
for the photometric error per frame. With these averages, we 
calculate the i iormalised Stokes parameters, QJI and Ujl with 
the formulae dZapatero Osorio etakl 120051 : iPatat & Romanieliol 
[2006): 



fo(n)-fE(n) 

Mrf) + fEiji) 
1 




Pq = QH = -^(Fm - F(45°)) 
Pu = UII = ^(F(22.5°) - F(67.5°)) 



4 5 

SNR/100 

Fig. 1: Polarimetric accuracy cr(PQ) (filled symbols) and cr{Pu) vs. 
signal-to-noise ratio per frame (SNR), in the / band (squares: 
P75, triangles: P77). The SNR is the averaged SNR over all 
extraordinary and ordinary fluxes for the given object. The 
(t{Pq) and cr(Pu) values are linked by a dashed line for a 
given target, for clarity (solid line for LHS102BC, dotted line 
for WD1344+106). The horizontal line indicates the theoretical 
limit for cr(Pi) of Patat & Romaniello (2006). The vertical thick 
line indicates the statistical l-cr uncertainty on the error expected 
for three measurements per retarder position. 



pl = V(e//)2 + (f///)2 



where foiif) and fEirf) are the ordinary and extraordinary aver- 
aged fluxes measured at the retarder angle rj. 

We propagate the photometric errors in order to obtain the 
errors on Qjl and Ujl (see Fig. [Til. Using the standard deviation 
distribution of multiple ^(77) measurements is a robust, model- 
independent way to estimate the polarimetric eiTors. In partic- 
ular, it is insensitive to eiTors in the CCD gain or in the noise 
model. However, it does not include flat-fielding errors, which 
in our case have a systematic effect for each target. Furthermore, 
the small number of images (generally 4x3) results in a large 
relative uncertainty (roughly 30%). The case of the polarisation 
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Table 3: Polarisation degrees and observing dates of standard stars: highly-polarised and (in italics) non- polarised. The Qll a nd t/// values are not 
corrected for the pleochroism (see text). The reference polarisation degrees and angles are from lFossati et alj ( |2007|) . 



Ohiert 


Date (\ IT) 


Filter 


o/i r%i 


[f/i r%i 


Pn, r%i 

^ Obs L '^i 


-* Foss 


r%i 


Onu (°) 




WD1615-154 


02/08/05 


R-Bessel 


-1-0.13 


+0.16 


0.14 






-25 




WU1j44+JUo 


A 1 /AO /AC 

01/08/05 


I-Bessel 


A A"7 1 A 1 Q 

-0.07 ± 0.13 


1 A 1 Q 1 All 

+U.13 ± 0.11 


0-0.199 






+34 ± 23 




WD03 10-688 


19/08/06 


1-Bessel 


+o.n 


+0.03 


0.17 






+5 




WD0310-688 


21/08/06 


I-Bessel 


+0.05 


+0.06 


0.08 






+24 




Hiltner 652 


03/08/05 


V-Bessel 


-F6.29 


+0.33 


6.30 


6.25 + 


0.01 


179.7 


179.18 ± 0.20 


Hiltner 652 


03/08/05 


R-Bessel 


+5.90 


-0.28 


5.91 


6.07 + 


0.02 


179.8 


179.39 ± 0.10 


NGC 2024 


03/08/05 


R-Bessel 


-0.22 


-9.63 


9.63 


9.62 + 


0.01 


45.5 


135.84 ± 0.02 


NGC 2024 


03/08/05 


1-Bessel 


-0.48 


-9.13 


9.14 


9.12 + 


0.04 


46.4 


136.26 + 0.15 


NGC 2024 


03/08/05 


V-Bessel 


+0.57 


-9.58 


9.60 


9.65 + 


0.06 


134.9 


135.47 ± 0.59 


Hiltner 652 


01/08/06 


I-Bessel 


+5.43 


-0.64 


5.47 


5.61 ± 


0.04' 


179.5 


179.18 + 0.11' 


Hiltner 652 


02/09/06 


I-Bessel 


+5.55 


-0.59 


5.52 


5.61 ± 


0.04' 


179.8 


179.18 + 0.11' 


NGC 2024 


03/09/06 


R-Bessel 


-0.11 


-9.58 


9.58 


9.62 + 


0.01 


45.9 


135.84 ± 0.02 


NGC2024 


10/08/06 


I-Bessel 


-0.59 


-9.16 


9.18 


9.12 + 


0.04 


46.0 


136.26 + 0.15 


BD-125133 


02/09/06 


I-Bessel 


+ 1.04 


-3.38 


3.54 


3.57 + 


0.09 


146.4 


143.99 ± 2.27 



' iFossati et al.l j2007i) only reports measurements using the PMOS mode. For the B and V filters for which they report measurements in both 
IPOL and PMOS modes, one finds difi'erences up to (0.10 + 0.02)% and (0.38 ± 0. 14) deg. 



degree is more complicated. The error propagation gives: 



p2 , pi 



Because it is the quadratic sum of two variables, that are 
roughly Gaussian distributed, its distribu tion is biased to- 
ward non-zero measure ments (see e.g. iFendt et al] Il996t 
iPatat & Romanielloll2006l and references therein). This prob- 
lem has been studied in a number of publications, and we 
adopt the correction of th e polarisation degree proposed by 
ISunmons & Stew^ (Il985l) . as well as their confidence inter- 
val determination — their Fig. 6, taking as normalised polarisation 

degree p = ^Pl/(r(PQ)^ + Pl/cr(Pu)^. 

When a significant polarisation degree is mea- 
sured, we ca lculate th e polaris ation angle, which 
is unbiased (IWardle & Kronbergi Il974 . We follow 
iLandi Degl'Innocenti et al.l (^07), whose expression takes 
care of the usual discontinuity around Pu = 0: 



- sign(f f/) arccos 



4 



Pl + Pl 



For significant polarisa tion degrees (Pi > 3o-(Pl), 
IWardle & Kronberg 1974), the bias on the polarisation an- 
gle becomes negligible and we simply propagate the errors: 



(r(Puf 



\Pq + PI/Pq\ 



FORSl sufi^ers from in strumental polarisation in the edges 
(iPatat & Romanielloll2006l) . but it is negligible at the centre of 
the chip. We repeatedly observed our targets at this central po- 
sition, so that our observations do not suffer from instrumen- 
tal polarisation. We measured the polarisation degree of stars 
in the field of view to check the instrumental polarisation, flat- 
fielding and interstellar polarisation. For those stars away from 



th e centre, we corrected the 
of iPatat & Romanielld ( [20& 
available in the z band. 



olarisation degree using the data 
in the / band. No such data are 



2.4. Polarimetric standard stars 

FORSl has proved to be extreinely s table in its polarimetric 
mode (M02; iPatat & Romaniellol2006h . The instrumental polar- 
isation is checked on a monthly basis with highly polarised stan- 
dard star observations, and more frequently with non-polarised 
standard stars. It appears that the polarisation angle needs to 
be corrected for instrumental polarisation by a constant (filter- 
dependent) amountfl while no correction in the polarisation de- 
gree is required. 

Our error determination procedure requires multiple images 
for each retard er angle, while we only ha ve one for the bright 
standard stars. iPatat & Romaniellq (l2006l) provide an estimate 
of the minimal uncertainties that can be obtained with the four- 
retarder-position procedure of 0.05% (see Fig.[T]i. Six of our po- 
larimetric measurements of polarised standard stars effectively 
fall within 1 cr of the published value, with one additional within 
2cr. 

The situation with non-polarised standard is not so clear: half 
the Q/I and U/I measurements are above 0.13% (in absolute 
value), while the other half satisfactorily lies within 0.07% of a 
null measurement. We find positive or negative measurements, 
as expected from random noise. The polarisation vector orienta- 
tion similarly does not exhibit any preferred direction. For one 
non-polarised standard star only, on August 1, 2005, we took the 
3-image-per-retarder-angle observing sequence also used for the 
science targets, and find a large polarimetric error of about 0.1% 
(see TableO. We could not find an obvious explanation for the 
large error Finally, the polarisation angle of NGC 2024 oscil- 
lates between 45 deg and 135 deg, which is due the discontinuity 
of the polarisation angle at Q/I - 0. This reflects a somewhat 
poor measurement of that parameter, while the polarisation de- 



^ The effect is called pleochroism jPatat & Romanieiioll2006h . For 
this chromatic zero angle, we used the values indicated in the FORS 1+2 
manual, section 4.6.2: -1.19deg inR, -2.89 deg in /, -1.64deg in z- 
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Table 4: Polarimetric results. The Qll and (/// values are not corrected for the pleochroism (see text). The last two columns give the 16% and 
84% percentiles of the P distributions. 



Object 


Date 


Filter 


an [%] 


UII [%] 


/>[%] 


^16% [%] 


^84% [%] 




DENIS J2200 


01/08/05 


I-Bessel 


+0.08 ± 0.09 


-0.01 ±0.08 








0.137 


3 


Kelu 1 AB 


02/08/05 
02/08/05 


I-Bessel 
z-Gunn 


+0.03 ± 0.05 
-0.02 ± 0.06 


+0.13 ±0.10 
+0.00 ± 0.07 


0.029 







0.201 
0.036 


2 
2 



2MASS J0036 


01/08/05 


R-Bessel 


-0.11 ±0.11 


+0.10 ± 


0.09 


0.065 


0.020 


0.222 


3 




Ul/Uo/U3 


I-Bessel 


— U.Uo ± U.Uj 


— U.U2 ± 


U.Uj 


n (Yin 
U.U/ / 


U.U4o 


n 1 (Yi 
U.iU/ 


i 




10/08/06 


I-Bessel 


-0.06 ±0.10 


-0.02 ± 


0.12 








0.112 


4 




11/08/06 


I-Bessel 


-0.03 ± 0.05 


-0.00 ± 


0.07 








0.054 


4 




average 


I-Bessel 


-0.06 ± 0.03 


-0.01 ± 


0.04 


0.053 


0.029 


0.079 




2MASS J2158 


01/08/05 


I-Bessel 


-0.01 ± 0.22 


-0.03 ± 


0.26 








0.057 


3 




02/08/05 


z-Gunn 


-0.08 ±0.11 


-0.19 ± 


0.19 








0.319 


3 


LHS 102BC 


01/08/05 


I-Bessel 


+0.38 ± 0.06 


-0.03 ± 


0.12 


0.376 


0.315 


0.435 


2 




02/08/05 


I-Bessel 


+0.23 ±0.15 


-0.08 ± 


0.22 


0.204 


0.057 


0.358 


2 




02/08/05 


z-Gunn 


-0.5 ±0.7 


-0.6 ± 


0.6 








1.209 


3 




78/07/06 




+0 30 + fl 33 


+0.22 + 


0.10 


0.321 


0.191 


0.499 


4 




average' 


I-Bessel 


+0.30 ± 0.05 


+0.04 ± 


0.07 


0.300 


0.246 


0.355 




2MASS J 1507 


01/08/05 


R-Bessel 


+0.10 ±0.15 


+0.04 ± 


0.27 








0.187 


3 




01/08/05 


I-Bessel 


+0.02 ± 0.05 


-0.00 ± 


0.11 








0.036 


3 


DENIS J2252 


02/08/05 


I-Bessel 


-0.07 ± 0.09 


-0.11 ± 


0.06 


0.116 


0.057 


0.181 


3 


2MASS J2244 


02/08/05 


z-Gunn 


+0.42 ± 0.84 


-0.26 ± 


0.20 








0.749 


2 








+0.32 ± 0.46 


-0.07 ± 


0.15 


02 





0.562 


2 




19/08/06 


I-Bessel 


+0.75 ± 1.46 


-0.52 ± 


1.08 








1.581 


3 








+0.04 ± 1.1 


-0.28 ± 


0.61 


0^ 





0.507 


3 




21/08/06 


I-Bessel 


-0.72 ± 0.88 


+0.82 ± 


0.87 








1.684 


3 








-0.54 ± 1.0 


+0.04 d 


: 1.1 


02 





0.967 


3 




22/08/06 


I-Bessel 


-0.39 ± 1.4 


+ 1.0± 


1.2 


02 





1.828 


3 




average'-^ 


I-Bessel 


-0.30 ±0.17 


+0.25 ± 


0.38 


0.3432 


0.157 


0.549 




6 Indi Bab 


01/08/05 


I-Bessel 


+0.03 ±0.13 


+0.05 ± 


0.13 


0.055 


0.000 


0.111 


3 



' Average based on the three /-Bessel measurements, the errors are propagated from the individual measurement uncertainties. 
2 PSF-photometry. For 2MASS J2244, only the three PSF-photometry results are considered in the average. 
^ Number (Ng) of exposures per retarder angle. 



gree fa lls within 0.05% of the published values of iFossati et al.l 
d2007h . 

All this leads us to believe that our minimal uncertainties 
probably lie around 0.1% rather than 0.05%, at least for some 
targets (see Fig.[T]i, which includes photonoise (we used very 
short integrations for the standard stars). In particular we do 
not believe our non-polarised standard stars measurements to 
be indicative of a systematic error on the polarisation degree or 
the polarisation angle. This 0.07% uncertainty is close or be- 
low the uncertainties on the measurements of our science tar- 
gets, and is probably a minor contributor compared to the tar- 
gets' photonoise. 

3. Results 

3.1. New polarimetric measurements 

We report our measurements in Table|4] All measurements of 
Q/I, U/I and Pi appear to be within 3cr of a null polarisa- 
tion, with a vast majority within 2cr. The only exception is the 
Q/I = -0.38+0.06%ofLHS102BConAugustl, 2005. Because 
this measurement has an uncertainty much smaller than usual 
for its SNR of 186 (see Fig.[rii, this uncertainty could be (statis- 
tically) underestimated, and the detection may be insignificant. 



More convincingly, the three Q/I measurements display a small 
dispersion of 0.06% for an average of {Q/I) - -0.31%. Two 
other targets, 2MASS J0036 and 2MASS J2244, have three I- 
band measurements, which are all compatible within 2 cr. This 
prompts us to also calculate the average of the Qs and Us. Again 
we use the standard deviation to determine the uncertainty. We 
also fail to detect significant polarisation on the averaged values 
of 2MASS J0036 and 2MASS J2244. Therefore we find only tar- 
get, LHS102BC, to show a significant polarisation. We measure 
a polarisation angle of -6.6 + 9.3 deg. 

Taking these results at face value leads to a fraction of po- 
larised ultra-cool dwarfs detected in the optical smaller in value 
than reported by M02 (about 50%) and ZO05 (29 + 9%). Our 
full sample of nine targets leads to a fraction of 11^^^% (1-cr), 
and smaller than 39% at 95% of confidence level (C.L.), but 
some of our measurements have large uncertainties. If we re- 
strict the sample to those targets with polarisation errors smaller 
than 0.2%, we find eight targets observed usually in the / band 
(but also once in the z band and twice in the R band), leading 
to an upper-limit of 28% (95% C.L.) of ultra-cool dwarfs with a 
polarisation degree larger than 0.5%, in agreement with previous 
works (M02,ZO05). 
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We note that we treat binaries as if they were single ob- 
jects. This is essentially true for e Indi BC and DENIS J2252, 
which have a small brightness ratio in the / band, so that 
only the warmer primary contributes to the flux. The cases of 
DENIS J2200 and of LHS102BC are diff^erent, being nearly- 
equal-brightness binaries. The contribution of each nearly equal- 
mass components to the observed linear polarisation of J2200 
and LHS102BC is unknown. Both components may contribute 
with a similar amount provided the fact they share very similar 
atmospheres, parallel rotation axis and nearly identical rotation 
rates. If only one of the objects of the pairs is polarised, our sen- 
sitivity to the polarised light of the latter would be decreased by 
a factor of two. 



3.2. Comparison with previous surveys 



Table 5: Comparison of our polarisation degrees with published mea- 
surements, in percent (1-cr interval or error). For each target 
and filter, the second line gives the 2-cr interval, correcting the 
published values following Simmons & Stewart ( 1985). 



As mentioned in Section l2n we preferentially selected targets 
with previous polarimetric measurements, in order to search 
for variability. Seven of our targets have previous / data, with 
2MASSJ0036 also having previous R data. We compare our 
measurements in Table|5] (see also Fig|2]i. In that table we use 
the best single-epoch measurement of ZOOS for 2MASS J0036 
inR. 

Our polarisation degrees along Q and U fall within 3cr of 
those of M02, also obtained with FORSl, with the exception 
of f q(LHS102BC). The observing set-up of the latter study is 
different, with 16retarder angles used, and only one image per 
angle. This allows a different error determination and a better 
removal of the instrumental polarisation and other systematics. 
The total integration time is similar (within » 50%) but the see- 
ing of our images is significantly worse. The PqS of LHS 102BC 
differ by 5cr. 

Compared to ZOOS we tend to find smaller polarisation de- 
grees (six times out of seven). Both observing set-ups are strictly 
identical, but the instruments differ The repetition of the pattern 
is unexpected. Variability could be proposed to explain a frac- 
tion of the observed changes, but cannot possibly explain all six 
measurements. Some formally-significant differences in Pq or 
f (/ are certainly due to an underestimated uncertainty when the 
errors on Pq and Pu widely diff'er in ZOOS (see also our Fig.lTJ. 
Eventually, only two measurements show robust discrepancies: 
Pu for 2MASS J1S07 and Pq for 2MASS J2244. 

Apart from intrinsic variability, differences in our measure- 
ments could be due to the different bandpasses, with FORS 1 / 
filter 7 2 n m bluer than CAFOS ' s. How ever, Sengupta & Krishan 
(1200 lb and lSengupta & Kwokl (l200Sb both predict an increasing 
polarisation degree with wavelength over 0.6-1 yum for a wide 
range of dust gr ain size (results are presented for 0.1-1/im). 
ISenguptal (l2003h presents similar conclusions over 0.7-0.9;um 
for particle sizes between 0.2 //m and 30yL(m. Under their hy- 
potheses, the difference in effective wavelength is not the right 
explanation to the different results. 

Among all three studies, although we have a sample of well- 
measured targets with multiple epochs relatively small, we do 
find three brown dwarfs with seemingly variable polarisation 
degree. This result needs to be confirmed with a more homo- 
geneous data set, but offers a promising prospect for the study of 
the dust cloud dynamics. 

3.3. Compiiation of existing data 

In order to increase the size of the sample, we collect mea- 
surements from the work of M02 and ZOOS. Again we re- 





r 


this paper 




iVlUZ 


Kelu 1 AB 


/ 


0-0.20 


0.80 ± 0.27 








0-0.30 


0.36-1.22 




2MASS J0036 


R 


0.02-0.22 


0.61+0.20 








0-0.32 


0.39-0.81 






I 


0.03-0.08' 


0.05 ± 0.06 


0.199 ±0.028 






0-0.11' 


0-0.13 


0.141-0.253 




J 
1 


u~u.uo 


L.JO ± U.JJ 














LHS 102BC 


I 


0.246-0.355' 




0.105 ±0.036 






0.191-0.409' 




0.032-0.171 


2MASS J 1507 


I 


0-0.04 


1.36 ± 0.30 








0-0.09 


1.01-1.69 




DENIS J2252 


1 


0.06-0.18 


0.62 ±0.16 








0-0.25 


0.38-0.84 




2MASS J2244 


I 


0.16-0.55'-2 


2.45 ± 0.47 








0-0.76 


1.90-2.98 








^,+FWHM 


^,+FWHM 


^,.±FWHM 






[nm] 


[nm] 


[nm] 


Filter 


R 


657 ± 150 


641 ± 158 




Filter 


I 


768 ± 138 


850 ± 150 


768 ± 138 



' Averaged value; ^ PSF-photometry; ^ Central wavelength and full- 
widths at half-maximum of the filter profiles, in nm. 



strict the sample to well-measured L-dwarf targets, with po- 
larimetric uncertainties smaller than 0.2%. For consistency we 
only consider /-band data. We count multiple non-detections 
(2MASS J0036 in our data and those of ZOOS) and multiple de- 
tections (LHS102BC in ours and M02's) as one measurement. 
In contrast, we count conflicting results (again 2MASS J0036, in 
M02) twice, which reflects the low significance of each individ- 
ual observations, or possible intrinsic variations. The fraction of 
polarised dwarfs which we derive really is the probability to de- 
tect polarised light using (various) non-ideal instruments, while 
the true fraction is necessarily higher (even assuming the mea- 
surements and uncertainties are correct). 

We have seven such objects in our sample, 17 in ZOOS and 
all eight targets of M02. We have one polarised target at the 3- 
cr level. M02 has five (including LHS102BC, already counted). 
ZOOS have five targets showing significant polarisation degrees. 
ZOOS do not correct for the intrinsic bias in the polarisation de- 
gree. However, their targets showing significant polarisation de- 
gree also display a significant polarisation in a single Stokes pa- 
rameter (in U/I), hence the detections are not due to that bias. 
The total is nine objects among 30. The fraction of polarised L 
dwarfs is therefore 30;';'*'% (1-cr errors). 

Restricting further the sample to targets with 0.1% uncer- 
tainties and defining polarised objects as those with polarisation 
degrees measured to be larger 0.3%, we obtain five targets in 
our sample, all eight of M02 and five in ZOOS including again 
2MASS J0036, for a total of 19. Only LHS102BC exhibit such a 
large polarisation degree (this paper). The statistical probability, 
given our errors, that it has a polarisation degree smaller than 
0.3% is about S0%, so it counts as 1/2 as polarised and 3/2 as 
unpolarised (M02). The upper- limit on the fraction of polarised 
ultra-cool dwarfs is now 16% (9S% confidence level). 
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4. Interpretation 

We have outlined in the section[T] an ambitious polarimetric 
study of brown dwarfs. As described above, given the small sig- 
nal and the diversity of observing set-up (with M02) and of in- 
struments (with ZO02), our results of polarimetric variability re- 
quire additional data for confirmation. We c an however attempt 
to com pare our results with the models of ISengupta & Kwokl 
(I2005h and with those of the atmospheric variability surveys. 
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4.1. Comparison with \Senaupta & Krishari \2001\) models, 
and correlation witii fiattening 

Despite our small sample size, we can compare the non- 
detection of polarisation for the five ultra-cool dwarfs with 
vsin/ measurements (see Table[T] and Fig|2]i. Four dwarfs 
have larger-than-average (projected) rotational velocities, with 
V sin / > 25 km/s. The inclination / is unknown, so how does 
the polarisation depend on it? At a given rotational velocity, 
the projected flattening of the dwarf decreases when sin / de- 
creases, i.e. when seen from the pole (in the extreme case of 
sin / = 0), and hence the polarisation degree decreases as well, 
linearly with sin /. However, at a given projected rotational ve- 
locity, the true rotational velocity would increase by (sin/)"', 
and the flattening would inc rease: Following the Eqs. 1 and 2 of 
ISenpupta & Krishanl (1200 lb . for most relevant eccentricities, i.e. 
for 0.04 < e < O.SQ, the difference of the equatorial and polar 
radii is roughly proportional to (sin /)"^, and therefore the polar- 
isation degree is nearly independent of sin /. 

Fo r a 1 .4-yum grain size. Figure 3 of ISengupta & Kwod 
(l2005h predicts a polarisation degree larger than 0.5% for a poly- 
tropic index n - I and rotation velocities of 25 km/s and higher, 
with a minimum at L4.5 and 0.5% polarisation degree. Our data 
strong ly exclude these pr e dictio ns, so that under the hypothe- 
sis of ISengupta & Kwokl ( |2005|) . either the typical grain size 
is larger flian 1 .4 jjm or the polytropic index is (unexpectedly) 
much larger. Alternatively, some of those hypotheses may not 
be applicable: a thicker cloud deck would increase the multi- 
scattering and reduce the net polarisation. We also stress that 
those predictions assume that the scatterers are spherically sym- 
metric (Simmons 1982). The polarisation due to scattering on 
non-spherical grains is barely understood and extremely difficult 
to predict. Non-spherical grains could induce much smaller po- 
larisation degrees. 

Finally, the polarisation predictions are based on atmo- 
spheric temperature-pressure profiles, as well as dust con- 
tent models, such as t hose of Marley (19 99, priv.com.), 
iBurrows et all ( l200lh and ICooper et all (l2003l) . Those models 
are found to reproduce reasonably well the spectral energy dis- 
tribution (SED) of L and T dwarfs, but the resulting polarisation 
prediction may be more sensitive to the input parameters than 
the SED. 

Because it is so difficult to predict the polarisation degree, 
it is equally diflicult to predict the fraction of brown dwarfs 
with a polarisation degree larger than some given value, which 
is what we measure. Instead we could compare the distribution 
of polarisation degrees, o r its upper-limits, with spectral type. 
ISen gupta & KwoM (1200 5^ predicts a larger polarisation degree 
between LO and LI and between L5 and L6. We do not detect po- 



The flattening /(«!), is related to the eccentricity e by the relation: 
= 1 - (1 - fY. Jupiter has a flattening of 0.065 corresponding to an 
observed eccentricity of 0.35. A 50-Mjup, 0.7-Rjup brown dwarf with the 
Jupiter rotational period of 9.9 h would have a predicted eccentricity of 
0.04. 




L1 L2 L3 L4 
Spectral type 

Fig. 2: Polarisation degrees as a function of spectral type for our sam- 
ple (filled circles), ZO05 (empty squares: published values; the 
error bars are corrected for the polarisation bias) and M02 
(em pty triangles), for the same targets. We superimpose the 
ISen gupta & Kwok (2005) predictions for various rotational ve- 
locities (starting from the bottom: 15, 25, 30 km/s). We also re- 
port V sin i measurements of the sources (in km/s) when avail- 
able. 



larisation for our targets of those spectral types: DENIS J2200, 
DENIS J2252, or 2MASS J1507, which has a projected rota- 
tional velocity of 27 km/s, a median value for brown dwarfs. 
This does not come in support of the hypotheses made by 
ISengupta & Kwod (l2005h . but clearly our statistics is very lim- 
ited. 



4.2. Comparison with atmospheric variabiiity studies 

Regarding the impact of cloud cover patchiness, our re- 
sults should be easier to compare with those of the spec- 
tro/photometric variability studies. Indeed, the linear polarisa- 
tion degree is expected to increase for larger physical scale of 
the surface heterogeneities, just as the variability signal due to 
the brown dwarf rotation or to the dynamical variations in the 
cloud cover. On the other hand, as noted by Sengupta (priv.com.) 
in ZO05, atmospheric variability may require thick clouds to be 
noticeable, and the clouds' thickness would result in multiple 
scattering and a reduced polarisation degree. 

On the observational side the matters are even more com- 
plicated. For the polarimetric studies we suffer from a theory 
limitation: we cannot relate a (detectable) polarisation degree to 
a physical size, while variability studies suffer from a lack of 
reproducibility: depending on the survey, the fraction of vari- 
able objects vary from 10% to 50%, as the different instru- 
ment sensitivity and stability, and classification methods, vary 
widely (Bailer-Jones 2005; Goldman 2005; Koen et al. 2005). 
Moreover, very few objects or even none (again depending on 
the definitions) have shown repeated signs of similar variabil- 
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ity, or a stable periodic variability, as expected for rotationally- 
driven variations. 

It is therefore futile to compare the fraction of polarised 
brown dwarfs and of variable ones. Among our targets several 
have monitored for photometric and/or spectroscopic variability: 

- Kelu 1 AB probab l y has the b est case for variability 
dClarke et al.ll2002l l2003t ICaball ero 2006), and we do not 
find it to be polarised, above the 0.3% level (Z O05 do). 

- 2MAS S J0036 has been detecte d as variable by iBerger et alj 
(l2005h in its radio emission, by lMaitil (l2007l) in t he R a nd / 
bands (not simultaneously), not by iGelino et al.l ( |2002|) and 
ICaballero et al. (2003), and we do not find it to be polarised 
above the 0. 13% level (M02 and ZO05, in R only, do). 

- 2MASS J1507 has shown some limited vari ability, over a 
few hours (Goldman et al, in prep. ),but not bv lKoenI (120031) . 
and we do not find it to be polarised above the 0.1% level 
(ZO05 do). 

- iMorales-Calderon et al. i (120061) found some evidence 
for variability in their Spitzer/IRAC monitoring of 
2MASSJ2244 at 4.5 yum. We do not detect polarisation 
above the 0.76% level, while ZOOS do. 

This compilation reveals a complicated and possibly con- 
tradictory picture. Apart from the possible and general over- 
confidence of authors in their results, we may conclude that the 
atmospheres of brown dwarfs could be very dynamic and vari- 
able, even at large physical scales. 



5. Conclusions 

We have obtained FORS 1 polarisation measurements in the R, I 
and/or z band for nine field ultra-cool dwarfs and binaries with 
spectral types of the primary from M9 to Tl. We detect a sig- 
nificant polarisation for LHS102BC, with Pl = (0.31 + 0.06)%. 
We failed to detect polarisation at the 3-cr level in the other 
targets, setting an upper-limit to the fraction of ultra-cool 
dwarfs with polarisation larger than 0.5% to 28% (95% con- 
fidence level). We tend to find smaller polarisation degrees 
than Zapatero Oso rio et al.l ( l2005l) . and larger uncertainties than 
iMenard et al.. (.2002 ). probably due in part to worse observing 
conditions and diff'erent observing strategy. Combining all three 
data sets, we set an upper-limit of 16% (95%C.L.) to the fraction 
of polarised ultra-cool dwarfs with polarisation degrees larger 
than 0.3%. 



Our small sarnple do es not reproduce the predictions of 
ISengupta & Kwod (l2005l) regarding the flattening-induced po- 
larisation of brown dwarfs of known rotational velocities. 
Possible explanations include the eff'ects of non-spheroidal 
grains, of multi-scattering, or of the underlying atmospheric 
models, be it the dust distribution or the grain size distribu- 
tion. Regarding the polarisation induced by heterogeneous cloud 
coverage, we fail to find a correlation between the polarisation 
measurements and the results of atmospheric variability studies. 
The discordant results of both polarimetric or variability stud- 
ies may be due to a complex dynamics of the ultra-cool dwarf 
atmospheres. 
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